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The synthesis and crystallography, magnetic susceptibility \, magnetization M, specific heat C p , 
in-plane electrical resistivity p and in-plane magnetic penetration depth measurements are reported 
for single crystals of ylPd2As2 (A — Ca, Sr, Ba) versus temperature T and magnetic field H. 
The crystals were grown using PdAs self-flux. CaPd2As2 and SrPd2As2 crystallize in a collapsed 
body-centered tetragonal ThCr2Si2-type structure (14/mmm), whereas BaPd2As2 crystallizes in the 
primitive tetragonal CeMg2Si2-type structure (P4/mmm), in agreement with literature data. The 
p(T) data exhibit metallic behavior for all three compounds. Bulk superconductivity is reported 
for CaPd2As2 and SrPd2As2 below T c — 1.27 and 0.92 K, respectively, whereas only a trace of 
superconductivity is found in BaPd2As2. No other phase transitions were observed. The x(T) and 
M(H) data reveal anisotropic diamagnetism in the normal state, with Xc > Xab for CaPd2As2 and 
BaPd2As2, and Xc < Xab for SrPd2As2- The normal and superconducting state data indicate that 
CaPd2As2 and SrPd2 As2 are conventional type-II nodeless s-wave electron-phonon superconductors. 
The electronic superconducting-state heat capacity data for CaPd2As2, which has an extremely 
sharp heat capacity jump at T c , are analyzed using our recent elaboration of the a-model of the 
BCS theory of superconductivity, which indicates that the s-wave gap in this compound is anisotropic 
in momentum space. 

PACS numbers: 74.70.Xa, 74.25.Bt, 74.25. Op, 74.25.N- 



I. INTRODUCTION 

The observation of high-temperature superconductiv- 
ity (SC) with transition temperatures T c < 38 K in doped 
122-typc iron arsenides such as in Ai_ 2; K a ;Fe2As2 (A = 
Ba, Sr, Ca and Eu) compounds has stimulated great 
interest in these materials. 1_15 The parent compounds 
AFe2As2 with A = Ca, Sr, and Ba are itinerant antifer- 
romagnetic (AF) semimetals that undergo a structural 
distortion from a tetragonal structure to an orthorhombic 
one on cooling below room temperature. The structural 
transition precedes a long-range AF itinerant spin den- 
sity wave (SDW) transition. Superconductivity in these 
compounds emerges upon suppression of the SDW tran- 
sition that can be achieved either by partial chemical 
substitutions at either the A, Fe or As sites or by appli- 
cation of external pressure. 11 The same phenomenology 
is found in the high-T c cuprates where the long-range AF 
order must be largely suppressed by doping prior to the 
emergence of superconductivity, but strong dynamic AF 
spin fluctuations must still be present. Thus, the iron ar- 
senides and the high-T c cuprates have the same generic 
phase diagram for the emergence of superconductivity, 
even though the cuprate parent compounds are AF insu- 
lators rather than SDW semimetals. n ~ 13 ' 16 ~ 18 

The partial substitutions at the Fe-site in AFe2As2 
by transition metals M in Ba(Fei_ a; Mi : )2As2 with M 
= Cr, Mn, Co, Ni, Cu, Ru, Rh, and Pd have re- 
cently been studied both theoretically and experimen- 
tally from the perspective of the degree of charge doping 
because of the changes in the magnetic and SC proper- 
ties caused by such substitutions. M2, 19-39 The substitu- 
tions for the Fe atoms by the 3d elements Co and Ni 



and the 4d elements Rh and Pd which have higher num- 
ber of 3d electrons than that of Fe are found to induce 
superconductivity. 6,21 ' 22,27 ' 28 However, no superconduc- 
tivity is induced by Mn or Cr substitutions having a lower 
number of 3d electrons. 23 ~ 26 On the other hand isoelec- 
tronic substitution of Fe by the Ad transition metal Ru is 
also found to induce the superconductivity 29,30 The com- 
mon feature in the SC materials is that the long-range 
SDW order in the parent compounds must be largely 
suppressed before SC occurs as noted above, where, in 
addition, the AF spin fluctuations are still strong. The 
latter fact suggests an unconventional magnetic mecha- 
nism for the SC in the FeAs-based materials. The extent 
and effects of electron and hole doping by such substitu- 
tions are still being debated. 12 ' 31 ' 37,39 ~ 41 

While the effect of partial substitutions of Fe by transi- 
tion metals have been investigated extensively, the phys- 
ical properties of the end-point compounds ^4M2As2 of 
the ^4(Fei_ a; M a ;)2As2 series are often not known in detail. 
Recently we reported crystallography and physical prop- 
erty studies of the end-point compounds SrCu2As2 and 
CaCui.7As2 for M = Cu which are sp-band metals. 40,41 
These compounds form in a collapsed tetragonal (cT) 
structure in which the formal oxidation state of As is 
As~ 2 = [As-As] _4 /2 due to the intcrlaycr As-As bonding 
in the cT phase. Thus with the formal oxidation states of 
Sr and Ca taken to be +2, the Cu atoms in these cT com- 
pounds are then Cu +1 with a nonmagnetic 3d 10 config- 
uration as observed 40 ' 41 and also as previously predicted 
from electronic structure calculations for SrCii2As2 by 
Singh. 42 Thus, extrapolating from high to low Cu con- 
centrations, the electron or hole doping effect resulting 
from substitution of M = Cu for Fe in ^4(Fei_ a; M ; ;;)2As2 
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FIG. 1: (Color online) (a) A unit cell of ThCr 2 Si 2 -type body- 
centered tetragonal crystal structure (14/mmm) of CaPd 2 As 2 
and SrPd 2 As 2 . (b) Two unit cells, stacked along the c-axis, 
of the CeMg 2 Si 2 -type primitive tetragonal crystal structure 
(P4/mmm) of BaPd 2 As 2 . 



and T c — 0.92(5) K for SrPd2As2 arc also reported to- 
gether with other superconducting-state properties de- 
rived from the above measurements as well as from mag- 
netic penetration depth measurements. We analyzed the 
superconducting-state electronic entropy and heat capac- 
ity data using our recent elaboration of the a-model of 
the BCS theory of superconductivity. 46 The data indi- 
cate that CaPd2As2 and SrPd2As2 are both conventional 
nodcless type-II electron-phonon superconductors, but 
with anisotropic s-wave gaps. Filamentary superconduc- 
tivity was also detected below 2.0 K in BaPd2As2- 

The experimental details are described in Sec. IF 
The crystallographic studies of CaPd2 As2 , SrPd2 As2 and 
BaPd2As2 are presented in Sec. Ill and the physical prop- 
erties in Sees. IV, V, and VI, respectively. Our mea- 
surements of the superconducting magnetic penetration 
depths versus temperature for CaPd2As2 and SrPd2As2 
are presented in Sec. VII. A summary of the experimen- 
tal results and analyses and the conclusions are given in 
Sec. VIII. 



compounds is ambiguous and interesting. 12 ' 39,40 

To provide insight into the T-x phase diagrams 
of A(Fei_ x M E ) 2 As2 systems with M = Pd, we 
report here our investigations of the crystallogra- 
phy and physical properties of single crystals of 
the three (Ca,Sr,Ba)Pd2As2 end-point compounds. 
CaPd2As2 and SrPd2As2 are reported to form in 
the body-centered tetragonal ThCr2Si2-type structure 
(space group IA/mmm) and BaPd 2 As 2 in the primi- 
tive tetragonal CcMg2Si 2 -typc structure (space group 
P4/mwm), 43 ' 44 as shown in Fig. 1 for CaPd2As2 and 
BaPd2As2- Both structures contain similar Pd square- 
lattices. However, while the Ca atoms in CaPd2As2 form 
a body-centered tetragonal sublattice, in the BaPd2As2 
structure they form a primitive tetragonal sublattice. 
Furthermore, the fourfold coordination of Pd by As is 
tctrahcdral in CaPd2As2 but is planar rectangular in 
BaPd 2 As2- The BaPd2As2 structure can be derived from 
that of CaPd2As2 by a [i, i, 0] translation of the central 
As-Ca-As layer. Superconductivity was recently reported 
below T c = 3.0 K in the similar compound SrPd 2 Gc 2 
with the ThCr2Si2-type structure. 45 

In this paper we report crystallographic investiga- 
tions of crushed single crystal powders and physical 
property measurements of single crystals of the three 
(Ca,Sr,Ba)Pd2As 2 compounds using magnetic suscepti- 
bility x, isothermal magnetization M, specific heat C p , 
and afr-plane electrical resistivity p measurements as a 
function of temperature T and applied magnetic field 
H. The published structures of these compounds are 
confirmed. The p(T), X (T) and C P {T) data demon- 
strate that all three compounds are sp-band metals. The 
x(T) data indicate anisotropic diamagnetic behavior in 
all three compounds with \ c > Xab for CaPd2As2 and 
BaPd2As2 and Xc < Xab for SrPd2As2. Bulk super- 
conducting transitions at T c = 1.27(3) K for CaPd2As2 



II. 



EXPERIMENTAL DETAILS 



Single crystals of APd2As2 (A = Ba, Ca, Sr) were 
grown using PdAs flux. The starting materials were 
high-purity elemental Ca (99.98%), Pd (99.998%) and 
As (99.99999%) from Alfa Aesar and Sr (99.95%) and Ba 
(99.99%) from Sigma Aldrich. The Ca, Sr or Ba and prc- 
reacted PdAs flux taken in a 1:4 molar ratio were placed 
in alumina crucibles and sealed inside evacuated fused 
silica tubes. The sealed samples were heated to 1100 °C 
at a rate of 60-80 °C/h and held there for 12 h. After 
cooling at a rate of 2.5 °C/h to 800 °C, shiny platc-likc 
crystals of typical size 2 x 1.5 x 0.4 mm 3 were separated by 
decanting the flux with a centrifuge at that temperature. 

The phase purity and the chemical composition of 
the CaPd2As2, SrPd2As2 and BaPd2As2 crystals were 
measured using a JEOL scanning electron microscope 
(SEM) equipped with an energy-dispersive x-ray analyzer 
(EDX). High- resolution SEM images demonstrated that 
the crystals were single-phase. The EDX analyses of two 
or three single crystals of each compound were consis- 
tent with the ideal 1:2:2 stoichiometry. The crystal 
structures of the samples were determined by Rietveld 
refinements of powder x-ray diffraction (XRD) data for 
crushed single crystals collected on a Rigaku Gcigcr- 
fiex x-ray diffractomcter using Cu K Q radiation. The 
Rietveld refinements were carried out using the FullProf 
software 47 which also confirmed the single-phase nature 
of the crystals. 

The magnetization measurements were performed us- 
ing a Quantum Design, Inc., superconducting quantum 
interference device magnetic properties measurement sys- 
tem (MPMS). We use Gaussian cgs units for the magnetic 
results and discussion, where the unit of magnetic field 
H is the Oc = G, but we also use the Tesla, with 1 T 
= 10 4 Oe, as a unit of convenience (the magnetic mo- 



merit output by the MPMS software is in Gaussian cgs 
units where 1 emu = 1 Gem 3 ). 11 In order to obtain the 
magnetic moment of the crystals in a particular field H 
we subtracted the pre-calibrated magnetic moment of the 
sample holder (quartz rod/plastic piece and GE varnish 
that were used to mount the sample for MPMS measure- 
ments) in the same H. Due to the small magnitudes of 
the magnetic moments of the small crystals and the un- 
certainty in the magnetic moment of sample holder, the 
accuracy of the M(H,T) and x{T) data of the samples 
reported here is of order 10%. 

The heat capacity measurements were carried out us- 
ing the relaxation method in a Quantum Design, Inc., 
physical properties measurement system (PPMS). Tem- 
peratures down to 0.45 K were obtained using a 3 Hc at- 
tachment to the PPMS. The a6-plane p(T) measure- 
ments were performed by the standard four-probe ac 
technique using the ac transport option of the PPMS. 
The electrical leads were 25 /xm diameter platinum wires 
attached to the crystals with EPO-TEK P1011 silver 
epoxy cured at 110 °C for one hour. The accuracy of 
p is of order 10% due to uncertainties in the geometric 
factor. 

The temperature variation of the magnetic penetration 
depth A ff was measured using a tunnel diode resonator 
(TDR) technique operating at about 15 MHz. The res- 
onator consists of an LC tank circuit with a single-layer 
coil of inductance L ~ 1 pK, a capacitor with capaci- 
tance C ~ 100 pF and a tunnel diode that is biased to 
the region of negative differential resistance, thus com- 
pensating the losses in the circuit. The ci rcui t, therefore, 
self-oscillates at a frequency / = (27iV LC)^ 1 . When 
a sample with magnetic susceptibility \ is inserted into 
the coil, the total inductance changes and the resonant 
frequency shifts accordingly by an amount which is pro- 
portional to x- The temperature dependence of the res- 
onance frequency shift A/(T) induced by changes in the 
sample's magnetic response is related to the magnetic 
susceptibility \ an d hence A e ff(T) by 48-50 



penetration depths arc the respective values in the ab- 
plane. 



III. 



CRYSTALLOGRAPHY 



The room temperature powder XRD patterns for 
crushed single crystals of CaPd2As2, SrPd2As2 and 
BaPd2As2 are shown in Figs. 2(a), 2(b) and 2(c), respec- 
tively. The weak unindexed peaks marked with stars for 
BaPd2As2 arise from the small amount of flux attached to 
the surface of the crystals prior to crushing them for the 
measurements. The Rietveld refinements confirmed the 
reported ThCr2Si2-type body-centered tetragonal struc- 
ture (space group IA/mmm) of CaPd2As2 and SrPd2As2 
and the CeMg2Si2-type primitive tetragonal structure 
(space group PA/mmm) of BaPd2As2. 43,44 

The refinement profiles for these structural models are 
shown in Fig. 2. During the final refinements, the thermal 
parameters B were set to B = and the occupancies of 
the respective atomic sites were fixed to unity, since there 
were no improvements in the quality of fits or significant 
changes in the lattice parameters or in the As c-axis posi- 
tion parameter z^ s upon making small changes in B and 
in the occupancies. The crystallographic and refinement 
parameters are listed in Tables I and II. For comparison 
the lattice parameters and za s values from the literature 
are also listed in Table I, where good agreement is found. 

The c/a ratio and the interlayer As-As distance 
c^As-As = (1 — 2zas)c for (Ca,Sr)Pd2As2 and 2zasC for 
BaPd2As2 are listed in Table I, where a "layer" is defined 
as an As-Pd-As slab in Fig. 1. The values of c?As-As 
for CaPd2As2 and SrPd2As2 arc close to the covalent 
single-bond distance of 2.38 A for As, 52 indicating that 
these two compounds have cT structures. 40 Therefore, as 
discussed in Sec. I, the formal oxidation state of Pd in 
these (Ca,Sr)Pd2As2 compounds is Pd +1 , which is the 
same formal oxidation state as for Cu in ThCi'2Si2-type 
SrCu 2 As 2 . 40 



A/(T) = -4tt X (T)G 

Acff(T) 



(1) 



-Gil 



L c s 



tanh 



Aoff(T) 



where G is a sample- and coil-dependent calibration con- 
stant and L e g is the effective sample dimension. 48 The 
value of G is determined experimentally by pulling the 
sample out of the coil in situ and measuring the to- 
tal associated frequency shift, A/o, so that AA e ff(T) = 
i cff Sf(T)/Af where AA e ff(T) = A cff (T) - A cff (T min ), 
$f{T) = f{T) — /(T m i n ) and T m i n is the minimum tem- 
perature of the measurement. 48 Another way to deter- 
mine G is by matching the temperature dependence of 
the skin depth, d(T), obtained from the resonator re- 
sponse in the normal state to the measured resistivity, 
p(T), by using the relation d{T) = (c/2Tr) v / p(T)/f. 51 
The ac magnetic field was applied along the c-axis of sin- 
gle crystals of CaPd2As2 and SrPd2As2, so the reported 



IV. PHYSICAL PROPERTIES OF CaPd 2 As 2 
CRYSTALS 

A. Electrical Resistivity 

The in-plane p of CaPd2As2 versus T is shown for crys- 
tal #1 at H = in Fig. 3(a) and for crystal #2 under 
various H in Fig. 3(b). A sharp zero-field superconduct- 
ing transition is observed in Fig. 3(b) at T c = 1.28(3) K 
that is suppressed by H = 0.14 T to below 0.45 K. For 
crystal #1 the superconducting onset is also at 1.8 K as 
shown in the inset of Fig. 3(a). A broad superconducting 
transition was seen in another crystal (#3) with an onset 
at 1.8 K and zero resistance at 1.3 K (not shown). Inter- 
estingly the p(T) of crystal #2 in Fig. 3(b) exhibits an 
upturn and peak before the resistivity drops due to super- 
conductivity which is suppressed by a field H > 0.08 T. 



TABLE I: Crystallographic and Rietveld refinement parameters obtained from powder XRD data for crushed APd 2 As 2 (A = 
Ba, Ca, Sr) crystals. Error bars for the last digit of a quantity are given in parentheses and literature references are given in 
square brackets. 





CaPd 2 As 2 


SrPd 2 As 2 


BaPd 2 As 2 


Structure 


ThCr 2 Si 2 -type 


ThCr 2 Si 2 -type 


CeMg 2 Si 2 -type 




body-centered tetragonal 


body-centered tetragonal 


primitive tetragonal 


Space group 


IA/mmm 


IA/mmm 


P4/mmm 


Formula units/unit cell (Z) 


2 


2 


1 


Lattice parameters 








a (A) 


4.2824(2) 


4.3759(1) 


4.3438(2) 




4.299(1) [43] 


4.383(1) [43] 


4.346(1) [43] 




4.283(1) [44] 


4.380(2) [44] 




c(A) 


10.0880(4) 


10.1671(3) 


5.7536(2) 




10.102(2) [43] 


10.179(2) [43] 


5.758(1) [43] 




10.093(1) [44] 


10.169(1) [44] 




c/a 


2.3557(6) 


2.3234(4) 


1.3246(4) 


Vm (A 3 ) 


185.01(1) 


194.69(1) 


108.56(1) 


As c-axis coordinate zas 


0.3763(3) 


0.3768(1) 


0.2705(4) 




0.3796(2) [43] 


0.3766(2) [43] 
0.3768(1) [44] 


0.2700(8) [43] 


As-As interlayer bond 


2.496(7) 


2.505(2) 


3.113(4) 


distance dAs-As (A) 








Theoretical density (g/cm 3 ) 


7.230 


7.682 


7.648 


Molar volume Vm (cm 3 /mol) 


55.71 


58.62 


65.38 


Refinement quality 








x 2 (%) 


11.8 


3.68 


8.09 


R P (%) 


12.7 


5.41 


7.39 


Rv/p (%) 


18.6 


7.23 


11.2 



TABLE II: Atomic coordinates obtained from the Rietveld 
refinements of powder XRD data for crushed y4Pd 2 As 2 (A — 
Ba, Ca, Sr) crystals. 



Atom 


Wyckoff 
symbol 


X 


V 


z 


CaPd 2 As 2 (IA/mmm) 










Ca 


2a 











Pd 


Id 





1/2 


1/4 


As 


4e 








0.3763(3) 


SrPd 2 As 2 (IA/mmm) 










Sr 


2a 











Pd 


Id 





1/2 


1/4 


As 


4e 








0.3768(1) 


BaPd 2 As 2 (PA/mmm) 










Ba 


la 











Pd 


2e 





1/2 


1/2 


As 


2h 


1/2 


1/2 


0.2705(4) 



The origin of this p(T) peak is not clear. For crystal #1 
the residual resistivity at 2 K before entering the super- 
conducting state is po = 43 pQ cm and the residual resis- 
tivity ratio is RRR = p(300 K)/p(2 K) « 3. 

In the following we fit the normal-state p(T) data 
by the Bloch-Griineisen (BG) model which describes 
the electrical resistivity pec (T) due to scattering of the 
charge carriers by longitudinal acoustic lattice vibrations 



in the absence of Umklapp scattering, given by 53,54 
T \5 j.e H /T x 5 



PBG(T/e R ) = 4ft 



e R 



(e x - 1)(1 -e~ x ) 



-dx, 

) 

(2a) 

where Or is the Dcbye temperature obtained from fitting 
resistivity measurements and ft is a material-dependent 
prefactor that is independent of T. The experimental 
p(T) data are fitted by 



p(T)=p + p(@ R )p BG (T/Q R ), 
where Eq. (2a) yields 

p(t/q r = i) = Ail [ — - 

Jo (e. x 



(2b) 



l)(l-e- x ) 



—da; (2c) 



0.946 464 ft. 



The integral in Eq. (2c) can be obtained in closed form 
but the expression is cumbersome. 

An excellent fit of the p(T) data by Eqs. (2) was ob- 
tained using the three independent fitting parameters po, 
p(@ R ) and 9 R . for 2 K < T < 300 K, as shown by the 
red curve in Fig. 3(a). The single parameter 0r. deter- 
mines the T dependence of the fit. While fitting the 
p(T) data we used the analytic Pade approximant func- 
tion of T/Or. in place of Eq. (2a) as developed in Ref. 54 
which accurately describes P&g(T / Q R ) and greatly sim- 
plifies least-squares fitting of experimental p(T) data by 
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FIG. 3: (Color online) (a) In-plane electrical resistivity p of 
a CaPd2As2 crystal #1 versus temperature T measured in 
applied magnetic field H — 0. The red curve is a fit by the 
Bloch-Griineisen model in the T range 2-300 K. Inset: Ex- 
panded plot of p(T) below 10 K to show the onset of supercon- 
ductivity at 2.0 K. (b) Expanded plot of p(T) of CaPd 2 As 2 
crystal #2 for 0.45 K < T < 2.5 K showing the supercon- 
ducting transition for different values of H , as listed, applied 
along the c-axis. 



FIG. 2: (Color online) Powder x-ray diffraction patterns of (a) 
CaPd2As2, (b) SrPd2As2 and (c) BaPd2As2 recorded at room 
temperature. The solid lines through the experimental points 
are the Rietveld refinement profiles calculated [(a) and (b)] for 
the ThCr2Si2-type body-centered tetragonal structure (space 
group I4/mmm), and (c) for the CeMg2Si2-type primitive 
tetragonal structure (space group P4/mmm). In (a), (b) and 
(c), the short vertical bars mark the fitted Bragg peak posi- 
tions. The lowermost curves represent the differences between 
the experimental and calculated intensities. The unindexed 
peaks marked with stars correspond to peaks from residual 
PdAs flux on the surface of the samples. 



the BG theory. The fitting parameters are summarized 
in Table III along with those of the other two compounds 
discussed below. The 1Z value is obtained from the fitted 
value of p(6r) using Eq. (2c). 



TABLE III: Parameters derived from Bloch-Griineisen fits to 
the resistivities p within the ab-plane of the listed single crys- 
tals obtained using Eqs. (2). Here po is the residual resistivity 
extrapolated to T = 0, 0r is the Debye temperature deter- 
mined from resistivity measurements, p(On) is the fitted value 
of p at T — 0r, and TL is obtained from the fitted value of 
p(Or) using Eq. (2c). The p and 1Z values do not take into 
account the systematic error of order f 0% arising from un- 
certainties in the geometric factor required to calculate p of a 
crystal from its measured resistance. The accuracy of 0r is 
not affected by this systematic error since this parameter is 
determined solely from the T dependence of p. 



Compound 



po 
(fiQ cm) 



Or 
(K) 



p(e R .) 

(fiQ cm) 



n 

(pD. cm) 



CaPd 2 As 2 
SrPd 2 As 2 
BaPd 2 As 2 



42.60(7) 
7.57(6) 
1.02(1) 



135(3) 
170(3) 
114(1) 



30.3(5) 
27.6(5) 
8.84(1) 



32.0 
29.2 
9.34 




50 100 150 200 250 300 

T(K) 

FIG. 4: (Color online) The heat capacity C p of a CaPd 2 As 2 
single crystal [crystal #2 of Fig. 3(b)] versus temperature T 
for 2.5 K < T < 300 K measured in zero magnetic field. The 
red solid curve is the fitted sum of the contributions from the 
Debye lattice heat capacity CV Debye (T) and predetermined 
electronic heat capacity y n T according to Eq. (5a). Inset: 
Cp/T vs. T for 0.45 K < T < 2.5 K. The dotted red vertical 
line indicates the T c . 
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FIG. 5: (Color online) (a) Heat capacity C p versus tempera- 
ture T of a CaPd 2 As 2 single crystal [crystal #2 of Fig. 3(b)] 
for 0.45 K < T < 2.0 K measured in different magnetic 
fields H applied along the c-axis. (b) C p /T versus T 2 for 
0.45 K < T < 3.2 K with H = and 0.15 T. The blue curve 
is a fit of the H = 0.15 T data for 0.45 K < T < 4.5 K 
together with the H = data for 1.5 K < T < 4.5 K by 
Eq. (3). 



B. Heat Capacity 

1. Overview of the Superconducting and Normal State Heat 
Capacity 

The C P (T) of a CaPd2As 2 crystal is shown in Fig. 4 
for 0.45 K < T < 300 K. As shown in the expanded plot 
of the low-T data in the inset of Fig. 4, a sharp jump 
is observed in C P (T) due to the occurrence of supercon- 
ductivity at T c = 1.27(3) K. The C P (T) data measured 
under various H are shown in Fig. 5(a). These data show 
that the T c decreases with increasing H and is suppressed 



to below 0.45 K by H = 0.15 T. 

Plots of Cp/T vs. T 2 for H = and 0.15 T are shown 
in Fig. 5(b). We analyzed the normal-state C P (T) data 
at 0.45 K < T < 4.5 K for H = 0.15 T and the data at 
1.5K<T<4.5KforiJ = together according to 



Cp(T) 
T 



7n + (3T 2 + 5T 4 



(3) 



where j n is the normal state Sommerfeld electronic 
heat capacity coefficient, (3 is the Debye T 3 -law lattice 
heat capacity coefficient and 8T A is a higher-order lat- 
tice contribution. A fit of the normal-state data by 



TABLE IV: The linear specific heat coefficients 7 n and the coefficients /? and S of the T and T 5 terms in the low-T heat 
capacity, respectively, and the density of states at the Fermi energy T>{E F ) for both spin directions for CaPd2As2, SrPd2As2 
and BaPd2As2 single crystals. The Debye temperatures 0d obtained at low T and for all T from heat capacity measurements 
and the Debye temperature 0r obtained from fitting electrical resistivity data, respectively, are also listed. 



Compound 


7n 

(mJ/molK 2 ) 


(mJ/molK 4 ) 


S 
(/iJ/molK 6 ) 


V(E F ) 
(states/eVf.u.) 


6 D (K) 
from low-T 


e D (K) 

from all T 


e R (K) 

from p(T) 


CaPd 2 As 2 
SrPd 2 As 2 
BaPd 2 As 2 


6.52(2) 
6.43(3) 
4.79(2) 


0.463(6) 
0.369(8) 
0.638(5) 


6.8(4) 
3.7(5) 
4.0(3) 


1.87(1) 
1.89(1) 
2.03(1) 


276(1) 
298(3) 
248(1) 


252(2) 
245(3) 
227(2) 


135(3) 
170(3) 
114(1) 



Eq. (3), shown as the blue curve in Fig. 5(b), gives 
7 n = 6.52(2) mJ/molK 2 , /3 = 0.463(6) mJ/molK 4 and 
S = 6.8(4) /iJ/molK 6 . We estimate the Debye tempera- 
ture Od from j3 using the relation 55 



Bi 



12n 4 Rn 
5/3 



1/3 



(4) 



where R is the molar gas constant and n = 5 is the 
number of atoms per formula unit (f.u.), yielding 0rj = 
276(1) K. 

The measured normal-state C P (T = 300 K) = 
123 J/molK in Fig. 4 is close to the expected classical 
Dulong-Pctit high-T limiting value Cy = 3nR = 15 R 
= 124.7 J/molK at constant volume V due to acoustic 
lattice vibrations. 55,56 Our normal-state C P (T) data for 
the T range 2.5 K < T < 300 K were fitted by the sum 
of the above electronic term j n T and the Debye model 
lattice heat capacity Cy Debye (T) per mole of atoms due 
to acoustic lattice vibrations according to 

C p (T) = ln T + nC p (T = e D )CvDcb y c(r/e D ), (5a) 

where 56 



CvDebye(7 1 /©D) — 9-R 



j. \ 3 rQn/T x^e* 



e D 



(e x - l) 2 



and 



dx 

(5b) 



CvDcbyc(T 1 /OD — 1) — 9R 



1 4 x 

x e 



o (e^l) 2 
1.381 730 R, 



dx (5c) 



The integral in Eq. (5c) can be obtained in closed form 
but the expression is cumbersome. When carrying out 
the least-squares fit of the experimental C P (T) data by 
Eqs. (5) we used the high-accuracy analytic Pade approx- 
imant for Cy Debye (T/Qd) in Eq. (5b) that we formulated 
in Rcf. 54 that greatly simplifies the fit. The fit was car- 
ried out using the fixed value j n = 6.52 mJ/molK 2 ob- 
tained above. Thus 0d as the only adjustable parameter. 
The fit yielded 9d = 252(2) K and is shown by the red 
curve in Fig. 4. This value of 0d is close to but slightly 
smaller than the value 0d = 276(1) K obtained above 
from analysis of the low-T C p data. Such differences are 
expected due to the T dependence of 0d- 54 ' 56 



We estimate the density of states at the Fermi level for 
both spin directions V{E F ) from 7 n using the relation 55 



7n = 7o (1 + A e _ p h), 



7o 



■V(E F ), (6) 



where 70 is the bare (LDA band-structure) Sommer- 
feld coefficient and A G - P h is the electron-phonon coupling 
constant. The latter parameter can be estimated from 
McMillan's theory 57 for the electron-phonon mechanism 
of superconductivity and is related to ©d and T c by 



A,-. 



- P h 



1.04 + ^*ln(e D /1.45T c ) 
(1 - 0.62^*) ln(e D /1.45T c ) - 1.04' 



(<) 



Here /x* is the repulsive screened Coulomb parameter 
having a value between 0.1 and 0.15 and is usually taken 
as n* = 0.13. With this value of fi* together with 
T c = 1.27 K and 9d = 276 K as determined above (Ta- 
ble IV), Eqs. (6) and (7) yield 



A,-, 



- P h 



0.474, 



7o 



4.42 



mj 
molK 2 ' 



(8) 



The relatively small value of A G _ p h implies weak-coupling 
superconductivity in CaPd2As2- The value of D(E-p) is 
obtained from the value of 70 and Eq. (6) to be T>(Ep) = 
1.87(1) states/eVf.u. for both spin directions. 

The T>(E-p) and po can be used to estimate the mean 
free path £ for conduction carrier scattering at low T. 
In a single-band quasi-free electron Fermi gas model the 
Fermi velocity (speed) vp is 55 



2*3 



Tr 2 h 



V F 



m* 2 V{. u 



■T>(E F ), 



(9) 



where m* is the effective mass of the current carriers, 
Vf.u. = Kcii/2 from Table I is the volume per formula 
unit and Ti is Planck's constant divided by 27r. Assuming 
to* = to„ where to p is the free electron mass, we obtain 



v F = 1.76 x 10 8 cm/s. 



(10) 



The mean free path t = v f t, where r is the mean free 
scattering time, is obtained from v F using 55 



= 3tt^ 



h 

e 2 po 



m vf 



(11) 



where h/e 2 = 4108 0. Using po = 34.6 /iSlcm (for crys- 
tal #2), m* = m c and the value of v F in Eq. (10), Eq. (11) 



gives 



1.52 nm, 



(12) 



which is only 3.6 in-plane lattice constants (Table I). 

The plasma angular frequency uj p of the conduction 
carriers can be estimated using the 3D quasi-frce-electron 
single-band relation 



1 / rn^vp 
3^2 



yielding 



o5 



"-',"; 



4ime 2 4(m*e) 2 (v F /hY 



m* 3ir 

Using m* = m c and vf from Eq. (10) gives 

-.16 



UJ 



p = 1.94 x 10 16 rad/s. 



(13) 



(14) 



(15) 



The superconducting London penetration depth in the 
clean limit at T ~ 0, Al(0), is a normal-state property 
given by 46,58 



Al(0) 



(16) 



where c is the speed of light in vacuum. The value of cj p 
in Eq. (15) gives 



A L (0) = 1.54 x 10~ 6 cm = 15.4 nm. 



(17) 



However, we will see in the following section that 
CaPd2As2 is in the dirty limit and not in the clean limit 
that was treated by BCS, where the actual penetration 
depth is much larger than Al(0). 

The normal-state parameters obtained from the above 
analyses of C P (T) for CaPd2As2 are summarized in Ta- 
bles IV and V together with those obtained below for the 
other two compounds discussed in this paper. 

The Debye temperatures Or in Table III obtained for 
the three compounds from analyses of the respective p(T) 
data are also listed in Table IV for comparison. The large 
discrepancies between the values of Od and Or for each 
compound indicate that the assumptions 54 of the BG the- 
ory arc violated in the APd2As 2 compounds. The BG 
theory ignores Umklapp carrier scattering. Furthermore, 
if the Fermi wave vector k-p is significantly smaller than 
the Debye wave vector kjy , one one might expect ©d and 
0R, to be significantly different. Furthermore, a mecha- 
nism such as electron-electron scattering in addition to 
the electron-phonon scattering assumed in the BG theory 
could contribute to the T dependence of p. 



phonon contribution f3T 3 + ST 5 from C P (T) is shown 
in Fig. 6(a). The large sharp jump AC C in C c at 
T c = 1.27(3) K indicates that the superconductivity in 
CaPd2As2 is a bulk effect, confirmed below. The ver- 
tical heat capacity jump indicated by the vertical dot- 
ted line in Fig. 6(b) yields AC C /T C = 7.4(2) mJ/molK 2 
and hence AC C = 9.4(3) mJ/molK. Using the normal 
state 7„ = 6.52(2) mJ/molK 2 from Table IV, one obtains 
AC c /-f n T c = 1.14(3) which is much smaller than the BCS 
value AC e /-f n T c = 1.43 in the weak-coupling limit. 46,58 
Because the bulk superconducting transition in C C (T) is 
very sharp, we infer that the reduction in AC c /j n T c from 
the BCS value is intrinsic. This reduction can happen if 
the superconducting gap (order parameter) is anisotropic 
in wave vector space, either from anisotropy in the gap 
on a single Fermi surface 46 or from multiple bands with 
distinct Fermi surfaces, each with a different isotropic 
or anisotropic gap, 59 as discussed in the context of the 
a-model below. 



The Value of a within the a-Model 

For simplicity, we analyze the zero-field superconduct- 
ing state thermodynamic data for CaPd2As2 within the 
framework of the so-called a-model of the BCS theory 
of superconductivity, 46,60,61 where a = A(Q)/k-BT c . In 
this model the normalized gap A(T) = A(T)/A(0) is the 
same as in the BCS theory which is calculated from the 
BCS gap equation 46,60 with a = cubcs where 



«BCS 



— TT C ^ E 



1.7639 



(18) 



and 7e ~ 0.5772 is Eulcr's constant. However, for calcu- 
lating the thermodynamic properties one uses a variable 
a parameter, which represents an inconsistency in the 
model, but which allows one to fit superconducting state 
thermodynamic data that deviate from the BCS predic- 
tions. The superconducting and normal state electronic 
entropies at T c are the same, as in the BCS model, so the 
superconducting transition is still second order with no 
latent heat. One can interpret the quantitatively deter- 
mined deviations from the BCS predictions in terms of 
other models and theories. For example, values a > «bcs 
can arise from the presence of strong electron-phonon 
coupling in contrast to the weak-coupling assumption in 
the BCS theory, 62 whereas a < cubcs can occur from gap 
anisotropy in momentum space 46 which we therefore as- 
sume is responsible for the reduced heat capacity jump 
in CaPd2As2 compared to the BCS theory prediction. 

The value of a is related to the heat capacity jump at 
T c by 46 



2. Superconducting State Properties 



The electronic contribution C C (T) to C p (T) of 
CaPd2As2 obtained at low T by subtracting the low-T 



AC C (T C 

J n T c 



AC (T C ) 






12 
7<(3) V«bcs 



BCS V^BCS 
2 



(19) 



1.426 



f— 

V «BCS 
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FIG. 6: (Color online) (a) Electronic contribution C c to the 
heat capacity versus temperature T of CaPd2As2. (b) C e /T 
versus T. The solid red curve in each figure is the theoretical 
prediction of the a-model in Eq. (21a) for a = 1.58. The 
BCS value is obcs ~ 1.764. The normal-state heat capac- 
ity is C cn = 7 n r with 7 n = 6.52 mJ/molK 2 . In (a), this 
normal-state contribution is the blue straight line and the 
extrapolation to low T is the dashed blue line. In (b), the 
normal-state behavior is C an /T = y n . 



where £(x) is the Ricmann zeta function. Inserting our 
normalized experimental heat capacity jump value into 
Eq. (19) gives 



a = 1.58(2), 



(20) 



which is significantly smaller than the BCS value of 1.764. 

The temperature dependence of the superconducting- 
state electronic heat capacity is calculated within the a- 



C cs (t) 



6a 3 



/(I - /) 




(21a) 



where t = T/T c , the normalized normal-state electron en- 
ergy is e — e/A(0), the normalized excited quasiparticle 
(electron and hole) energy is 



E = 



E 

A(0) 



= V^ + A2, 



(21b) 



and the Fermi-Dirac distribution function in the dimen- 
sionlcss variables is (with E-p = 0) 



f = f(a,E,t) 

E 



1 



with 



knT 



e aE/t + I 

uE 

~T~' 



(21c) 



The t-dependent A(t) and dA(t)/dt 2 are calculated as 
described in Ref. 46. 

The solid red curves in Figs. 6(a) and 6(b) are the theo- 
retical predictions for C es (T) and C es (T)/T, respectively, 
calculated from Eq. (21a) using a = 1.58 and are seen to 
be in good agreement with the data. The agreement of 
the lowest-T data with the theory indicates that there is 
no residual electronic specific heat for T — > 0, which in 
turn indicates that the entire sample is superconducting 
with a single nodcless s-wave gap and with a single T c . 



Thermodynamic Critical Field H c 

The experimental thermodynamic critical field H c ver- 
sus T of a superconductor can be estimated using the 
zero- field C e (T) data via the electronic entropy differ- 
ence between the normal (S cn ) and superconducting (S cs ) 
states per unit volume at H = according to the 
Clausius-Clapeyron-likc relation 58,63 



S cn {T) - S es (T) 



1 dHl{T) 

8tt dT ' 



(22a) 



r-T', 



which, with S C {T) = f* [C c {T")/T")]dT", yields 

/.T c 

H 2 C {T) = 8tt / [S en (T') - S cs (T')]dT>. (22b) 

Jt 

In cgs units, H c is expressed in units of Oe, where 
1 Oe 2 = 1 erg/cm 3 , so S en (T') = j n vT' is in units of 
erg/cm 3 K. The Sommerfeld coefficient 7 n y in units of 
erg/cm 3 K 2 is calculated from the above Sommerfeld co- 
efficient 7„ in units of mJ/mol K 2 = 10 4 erg/mol K 2 ac- 
cording to 



7nV 



7n 

Vm 1 



(23) 
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FIG. 7: (Color online) (a) Thermodynamic critical field H c 
versus temperature T obtained for CaPd2As2 from the ex- 
perimentally derived electronic heat capacity C C (T) data us- 
ing Eqs. (22) (open circles) and the theoretical prediction 
of the a-model 46 ' 61 using T c = 1.27 K, a = 1.58 and 
7„ = 6.52 mJ/molK 2 in Eqs. (24) and (25) (red curve), 
(b) Upper critical magnetic field H C 2(T) of CaPd2As2 de- 
termined from the electrical resistivity p(T, H) (black open 
squares) and heat capacity C P (T,H) (blue open circles) data 
in Figs. 3(b) and 5(a), respectively. The red curve is the 
prediction for H c2 {T) by the WHH theory 65 in Eqs. (30) 
for dH c2 (T)/dT\ T = Tc = -0.18 T/K, a M = 0.13, and A so = 0. 



where Vm is the molar volume in units of cm 3 /mol. For 
CaPd2As2, our values j n = 6.52 mJ/molK 2 and Vm = 
55.7 cm 3 /mol from Table I give 7„v = 1170 erg/cm 3 K 2 . 
The experimental H C (T) data obtained from Eqs. (22) 
are plotted as open circles in Fig. 7(a). 

In the a-model, the thermodynamic critical field at 



T = is given by 



4(i 



flc(0) 



(7n V T?) 



1/2 



6 

— a 

■n 



1.382 a. 



(24) 



Using our values of 7„v, a = 1.58, and T c = 1.27 K gives 
H c (0) = 9.5 mT. The T dependence of H c is calculated 
from 46 






4tt 2 
3a^ j t 

47T 2 ^ 

3a 2 



&en\t ) &cs\t ) 



dt' 



t' 



S es (t') 



l-awT c 



dt' 



(25a) 



where the 

S'en(£')/7nvT , c 



normal-state electronic entropy is 
= £', and the superconducting state 



entropy is obtained from 
S es (t) 



ln\T c 



6a 2 



f(a,E,t)[E+-jr) 



(25b) 



The resulting H C (T) calculated for the above values of 
-ff c (0), a, 7 n v and T c is shown by the red curve in 
Fig. 7(a). Good agreement is observed between the tem- 
perature dependence of the data and the theoretical pre- 
diction of the a-model, although the calculated zero- 
temperature value appears to be a bit low. 



Upper Critical Field H C 2 

The p(T) data in Fig. 3(b) and the C P {T) data in 
Fig. 5(a) for CaPd 2 As2 in fields aligned along the c axis 
yield the T-dependence of the upper critical field H C 2 
shown in Fig. 7(b). One sees that H C 2(T — s- 0) is more 
than an order of magnitude larger than H c (0) = 9.5 mT 
determined above, indicating that CaPd2As2 is a type-II 
superconductor. Although the demagnetization factor of 
the plate- like crystal is large for H \\ c, this has no in- 
fluence on the present discussion because we only discuss 
the high-field behavior at the boundary with the normal 
state, where the magnetization is small. 

For a one-band typc-II BCS superconductor the orbital 
critical field H° 2 rh at T = is given by 64 ' 65 



tf c ° 2 rb 



(0) 



-AT, 



dH c2 (T) 



dT 



(26) 



T=T C 



where A = 0.73 and 0.69 in the clean and dirty lim- 
its, respectively. From Fig. 7(b), dH c2 (T) / dT\ T=Tc = 
-0.22(6) T/K as determined from the p(T) data for 
0.8 < T/T c < 1.0. Thus for CaPd 2 As 2 we estimate 
#c° rb (0) = 0.20(6) T in the clean limit and H% lb (0) = 
0.19(5) T in the dirty limit. 

The Pauli-limiting upper critical field at T = 0, iJp(0), 
is the maximum possible upper critical field and is the 
field at which the magnetic field energy of the current 
carriers is equal to the superconducting condensation 



11 



energy. 66,67 Within the a-modcl, this is given for spec- 
troscopic splitting factor g = 2 by 

^M = « w L2473 (^_) , (27) 

where /xb is the Bohr magneton and obcs is given in 
Eq. (18), yielding 



H P (0)[T] = 1.86T C [K; 



( a 



V ckbcs 



(28) 



Taking T c = 1.27 K, a = 1.58 from Eq. (20) and 
a BC s = 1-7639 from Eq. (18) gives H P (0) = 2.12 T. 
Since the measured H c2 (0) ~ 0.15 T - 0.07ff P (0) [sec 
Fig. 7(b)], the effects of Pauli limiting on H c2 should be 
minimal. The Maki parameter qm expresses the relative 
magnitudes of the orbital and Pauli limiting H c2 values 
as 68 



«M 



_ m gc° 2 rb (o) 



(29) 



which gives «m = 0.13 for CaPd2As2. 

To include the influence of Pauli limiting and spin- 
orbit scattering of quasiparticles on H c2 , we analyzed the 
H C 2(T) data within the Werthamer, Helfand and Hohen- 
berg (WHH) theory for a one-band type-II dirty-limit 
superconductor which calculates H c2 in terms of the or- 
bital, spin-orbit scattering and Pauli spin paramagnetism 
contributions in dimensionless variables as 65 



In 1 V [ 1 

t ^ I I2I/+1I 



|2i/+l| 



(aMh/tf 



(30a) 



|2i/ + l| + (A + Aso)/t 



where t = T/T c , A so is the spin-orbit scattering parame- 
ter and 



_4_\ H c2 (T)/T c 



tt 2 J dH c2 (T)/dT\ T=Tc 



(30b) 



WHH state that the applicability of their theory to 
a specific superconductor can be tested by comparing 
the value of om calculated from Eq. (29) with the value 
obtained from their alternative expression 65 



a M = -0.528 



dH c2 {T) 



dT 



T=T a 



(31) 



where the derivative is in units of T/K. Using our value 



dH c2 {T) 



dT 



T=T a 



-0.22(6) 



T 
K 



(32) 



obtained from Fig. 7(b), Eq. (31) gives a M = 0.12(3), 
which is indeed in agreement with the above estimate 
% = 0.13, thus indicating that the WHH theory is 
appropriate for fitting our H c2 (T) data. Reasonable 



agreement of H c2 (T) with the WHH theory prediction 
for cum = 0.13 was obtained from Eq. (30) with A so = 
and dH c2 (T)/dT\ T = Tc = -0.18 T/K as shown in 
Fig. 7(b), where H c2 (0) = 157 mT. The same value 
of dH c2 (T) / dT\T=T c is obtained from the dirty-limit 
relation 69 dH c2 {T)/dT\ T=Tc = 4.48 x 10 4 7„ v po, where 
p is in units of f2 cm, indicating that CaPd2As2 is in the 
dirty limit as further documented below. 



Ginzburg- Landau Parameter 

The Ginzburg-Landau parameter kql can be esti- 
mated from the relation 58 



H, 



KGL 



c:2 



V2H C 



(33) 



Using H c2 = 157 mT and H c = 9.5 mT gives kql = 
11.7 3> l/v2, characterizing CaPd2As2 as a type-II su- 
perconductor. An estimate of the lower critical field H c \ 
is then obtained from 58 



-Hci = H c 



hi k gl 
V2 



KGL 



(34) 



which for F c (0) = 9.5 mT and k gl = 11.7 gives H cl (0) = 
1.4 mT. 

The Ginzburg-Landau coherence length at T = 0, £(0), 
can be estimated from 58,63 



Bca(O) 



$n 



2<(0) 5 



(35) 



where $o = 2.07 x 10 7 G cm 2 is the flux quantum. Using 
H c2 = 157 mT this gives 



£(0) = 45.8 nm. 



(36) 



This £(0) is much larger than the mean free path £ = 
1.52 nm in Eq. (12) which indicates that CaPd2As2 is in 
the dirty limit. 

In the dirty limit the penetration depth at T — > in 
the notation of Tinkham is given by 58 



1 + ^ (dirty limit). (37) 



The T dependence of £ is given by 58 

£(T) it H c (0) A L (0) 



Co 2\/3ffc(T)A e ff(T) : 



(38) 



yielding the relationship between the zero-temperature 
Ginzburg-Landau £(0) and the BCS £o as 



m * A L (0) 



Combining Eqs. (37) and (39) gives 



(39) 



io 



2«/3(l + & 



(dirty limit). (40) 
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FIG. 8: (Color online) Zero-field-cooled magnetic susceptibil- 
ity x of a CaPd2 As2 single crystal versus temperature T mea- 
sured in a magnetic field H — 3.0 T applied along the c-axis 
(Xc H || c) and in the ab-plane (Xab, H _L c). The filled stars 
represent the intrinsic \ obtained from fitting M(H) isotherm 
data in the Appendix by Eq. (A. 2a) which are more accurate 
that the temperature-scanned data at fixed field. The lines 
joining the stars are guides to the eye. 



Substituting the above values for £(0) in Eq. (36) and £ 
in Eq. (12) into Eq. (40) and solving for £ gives £o = 
1690 nm. Then using Al(0) = 15.4 nm, £o = 1690 nm 
and £ = 1.52 nm. Eq. (37) gives A e g(0) = 510 nm in the 
dirty limit, which is of the same order as the value of 
210(60) nm obtained from the penetration depth mea- 
surements in Sec. VII below. 

The BCS coherence length is related to the Fermi ve- 
locity within the a-model by 46 



io 



frv-p 
ttA(O) 



1 



hvp 

knTc 



(41) 



This allows an additional estimate of the Fermi veloc- 
ity which for a = 1.58 and £o = 1690 nm yields vp = 
1.40 x 10 8 cm/s. This value of wp is close to the value 
v F = 1.76 x 10 8 cm/s in Eq. (10) obtained from V(E F ), 
indicating the self-consistency of our modeling. 

A summary of the measured and derived superconduct- 
ing state parameters for CaPd2As2 is given in Table V. 



C. Magnetization and Magnetic Susceptibility 

Zero-field-cooled (ZFC) x( T ) = M{T)/H data for a 
CaPd2As2 crystal versus T in a magnetic field H = 3.0 T 
applied along the c-axis (x c , H \\ c) and in the aft- 
plane (Xab, H _L c) are shown in Fig. 8. The data 
for both directions of H are strongly diamagnetic and 
nearly independent of T except for the Curie-like up- 
turns at low-T which were found from analysis of M(H) 



TABLE V: Measured and derived superconducting and rel- 
evant normal state parameters for CaPd2As2. T c : bulk su- 
perconducting transition temperature; y n : observed Sommer- 
feld coefficient of the linear term in the low-T normal-state 
heat capacity; A e - P h: electron- phonon coupling constant; 70: 
bare Sommerfeld coefficient; AC C : heat capacity jump at T c ; 
a — A(0)/kBT c ; A: superconducting order parameter; qm: 
Maki parameter; H c , Hp, H c i, H®2 h , H c %: thermodynamic, 
Paul limiting upper critical, lower critical, orbital upper crit- 
ical, and fitted upper critical magnetic fields, respectively; 
kgl: Ginzburg-Landau parameter; £: Ginzburg-Landau co- 
herence length; £0: BCS superconducting coherence length; I: 
electronic mean- free path at low T; ui p : angular plasma fre- 
quency; Al: London penetration depth; A c fi: magnetic pen- 
etration depth. The value of A°g s (0) is determined from the 
magnetic penetration depth measurements. 



CaPd2As2 property 



value 



Tc(K) 


1.27(3) 


7n (mJ/molK 2 ) 


6.52(2) 


A c — ph 


0.474 


70 (mJ/molK 2 ) 


4.42 


AC C (mJ/molK) 


9.4(3) 


ACo/7nT c 


1.14(3) 


a (from AC* c / 7n T c ) 


1.58(2) 


A(0)/fc B (K) (observed) 


2.02(14) 


am 


0.13 


H C (T = 0) (mT) 


9.5 


H F (T) 


2.12 


H e i(T = 0) (mT) 


1.4 


H° 2 Th {T = 0) (dirty limit) (T) 


0.19(5) 


# c2 (T = 0) (mT) 


157 


KGL 


11.7 


€(T = 0) (nm) 


45.8 


Co (nm) 


1690 


£ (m* = m e ) (nm) 


1.52 


uj p (m* — m c ) (10 16 rad/s) 


1.94 


A L alc (0) (clean limit) (nm) 


15.4 


A^ lc (0) (dirty limit) (nm) 


510 


A^ s (0) (nm) 


210(60) 



isotherms in the Appendix to be due to the presence of 
a small amount of saturable paramagnetic (PM) impu- 
rities. The intrinsic x(^) values at several temperatures 
obtained from the analysis of the M(H) isotherms are 
shown by filled stars in Fig. 8, which are more accurate 
than the x( T ) = M(T)/H data. The x is anisotropic 
with Xc > Xab over the entire T range. 

The x anisotropy in CaPd2As2 is different from the 
anisotropy usually observed in doped and undoped FeAs- 
bascd ThCr2Si2-structure compounds, where Xab > Xc 
such as in BaFe2As 2 . 10,11 The powder and temperature 
average of the intrinsic x obtained from fitting the M{H) 
isotherm data in the Appendix over the T range 20 to 
300 K is ( X ) = [2{ X ab) + (Xc)]/3 = -1.5 x lO" 4 cm 3 /mol. 

The diamagnetic susceptibilities in Fig. 8 seem to be 
rather large. For the purpose of comparing the dia- 
magnetic susceptibilities of different materials, a com- 
parison of the thermodynamic dimcnsionless magnetic 
susceptibility per unit volume xv (volume susccptibil- 
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ity) is most appropriate. Using the molar volume Vm = 
55.7 cm 3 /mol for CaPd 2 As 2 in Table I and Xab(100 K) = 
— 1.6x 10~ 4 cm 3 /mol from Fig. 8 (black star), one obtains 
Xv at = Xab/Vu in the ab plane as 

Xvah(100 K) = -0.29 x 1(T 5 (CaPd 2 As 2 ). (42) 

This value can be compared with the respective values 
for elemental Bi and C (graphite) that are well-known 
for their exceptionally strong diamagnetism. Bi has 
a rhombohedral crystal structure with a mass density 
p m = 9.8 g/cm 3 and a gram susceptibility in the hexag- 
onal ab plane x g a&(100 K) = -1.9 x 1CT 6 cm 3 /g. 70 The 
volume susceptibility xv — PmXg is then 



Xvab(100 K) = -1.9 x 1CT 5 (Bi). 



(43) 



On the other hand, using p m = 2.27 g/cm 3 for hexag- 
onal highly-oriented pyrolytic graphite (HOPG) and 
X gc (100 K) = -2.3 x 10 -5 cm 3 /g, 71 one obtains 



Xvc(100 K) = -5.2 x 10~ 5 (C, HOPG). 



(44) 



These values for Bi and graphite are about 6.5 and 18 
times more diamagnetic than the value for CaPd 2 As 2 in 
Eq. (42), respectively. However, these three values are 
all much less diamagnetic than the value xv = — 1/4-7T ~ 
—0.0796 for the diamagnetic susceptibility of a supercon- 
ductor at low fields with zero demagnetization factor due 
to complete exclusion of the magnetic induction from the 
interior (except within a magnetic field penetration depth 
of the surface). 

The different contributions to the intrinsic x are 



X = Xcoro + XVV + XL + XP> 



(45) 



where the first three terms are orbital susceptibilities and 
the last term is the Pauli spin susceptibility of the con- 
duction carriers. Xcoro is the isotropic diamagnetic sus- 
ceptibility of localized core electrons, xvv is the generally 
anisotropic paramagnetic Van Vlcck susceptibility, and 
Xl is the generally isotropic Landau diamagnetic suscep- 
tibility of the conduction carriers. 

The Xcoro is estimated using atomic diamagnetic 
susceptibilities 72 which gives Xcoro = —1-78 x 10~ 4 
cm 3 /mol. The xp is related to V(E F ) by 11 ' 73 



Xp - 9 -p 2 B V{E ¥ ). 



(46) 



Then using g = 2 and V(E F ) = 1.87 statcs/cVf.u. for 
both spin directions, from Table IV we obtain xp = 6.0 x 
10~ 5 cm 3 /mol. The xl is related to xp by 73,74 



XL 



1 / m c \ 
3 \m*J 



Xp, 



(47) 



where m* is the effective mass and m is the electron 
mass. Assuming m* = m e wc obtain xl ~ —2.0 x 
10~ 5 cm 3 /mol from the above value of xp- Then xvv 
is obtained by subtracting these three contributions from 



TABLE VI: Estimated contributions to the intrinsic angle- 
and temperature-averaged magnetic susceptibilities (x) of 
APd2As2 (A = Ba, Ca, Sr) crystals. Here xp is the Pauli 
spin susceptibility of the conduction carriers, and the orbital 
susceptibility contributions are the diamagnetism Xcore of the 
atomic electron cores, the Landau diamagnetism xl of the 
conduction carriers and the Van Vleck paramagnetism Xvv- 
All susceptibilities are in units of 10~ 5 cm 3 /mol. Possible rea- 
sons for the unphysical negative value of (xvv) for CaPd2As2 
are discussed in the text. 



Compound 


(x) 


Xcorc 


Xp 


XL 


(xw) 


CaPd 2 As2 
SrPd 2 As 2 
BaPd 2 As 2 


-15.4 
-3.9 
-12.4 


-17.8 
-19.3 
-21.6 


6.0 

6.0 
6.6 


-2.0 
-2.0 
-2.2 


-1.6 

11.4 
4.8 



the measured x according to Eq. (45). The four contribu- 
tions to the intrinsic x are summarized in Table VI, along 
with corresponding values for SrPd 2 As 2 and BaPd 2 As 2 
determined below. 

It is seen in Table VI that the inferred value of xvv 
for CaPd 2 As 2 is negative, which is unphysical. The rea- 
son for this error is not clear. The most likely source of 
the negative xvv value is a small error in correcting the 
total measured magnetic moment for the sample holder 
contribution, which was up to 40% of the total measured 
moment. Thus the uncertainty in the measured x values 
is of order 10% as discussed in Sec. II. In particular, the 
negative xvv is about 8% of the measured moment, and 
an error of only 4% in the sample holder correction could 
cause the derived xvv to be negative. 



V. PHYSICAL PROPERTIES OF SrPd 2 As 2 
CRYSTALS 

A. Electrical Resistivity 

The in-plane p(T) data for SrPd 2 As 2 measured on 
two different crystals #1 and #2 at different H are 
shown in Fig. 9. Metallic behavior is evident from the 
T dependence of p in Fig. 9(a). The expanded plot of 
p(T) in Fig. 9(b) reveals a superconducting transition at 
T c = 1.5(1) K. The T c is suppressed with increasing H, 
as shown. The data in Fig. 9(b) are noisy due to the small 
sample size and the small voltage signal arising from the 
small magnitude of the resistivity at low T. 

The two crystals were found to have different residual 
resistivities and RRR values. For crystal #1, /Co(L8 K) = 
7.5 pfl cm and RRR = p(300K/p(l.8K) w 8, whereas 
for crystal #2, po — 4.9 pft cm and RRR sa 4. Since 
crystal #1 has the higher RRR, we analyzed the p(T) 
data of this crystal in Fig. 9(a) using the BG model. A fit 
of the p(T) data in Fig. 9(a) by Eqs. (2) for 1.8 K < T < 
300 K gives p = 7.57(6) fj.fl cm, p(0 R ) = 27.6(5) fj.fl cm 
and Or = 170(3) K, where wc used our analytic Pade 
approximant in place of the integral in Eq. (2a) . The 
good fit obtained is shown by the red curve in Fig. 9(a). 
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FIG. 9: (Color online) (a) In-plane electrical resistivity p of 
a SrPd2As2 crystal (crystal #1) versus temperature T mea- 
sured in applied magnetic field H = 0. The red curve is a fit 
by the Bloch-Griineisen model in Eqs. (2). (b) Expanded 
plot of the low-T p(T) data for SrPd 2 As 2 crystal #2 for 
0.45 K < T < 4 K showing the superconducting transition 
for different H applied along the c axis. The noise in the data 
is due to the small size of the crystal and the small magnitude 
of p at low temperatures. 



The value of 1Z obtained from the value of p(6r) using 
Eq. (2c) is 1Z = 29.2 pil cm. The parameters obtained 
from the fit arc summarized in Tabic III. 



B. Heat Capacity 

The C P (T) data for SrPd 2 As 2 arc shown in Fig. 10. 
The C P (T = 300 K) « 124 J/molK is close to the classi- 
cal Dulong-Pctit high-T limit C v = 5R= 124.7 J/molK. 
In order to correlate low-T C P {T) data with p(T) data 




0.0 0.5 1.0 1.5 2.0 
T(K) 




100 150 200 250 300 
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FIG. 10: (Color online) The heat capacity C p of a SrPd2As2 
single crystal versus temperature T measured in zero mag- 
netic field H . The red curve is the fitted sum of the con- 
tributions from the Debye lattice heat capacity CvDebye(T') 
and predetermined electronic heat capacity *y n T according 
to Eq. (5a). Inset: Expanded plot of C p /T versus T for 
0.45 K < T < 2.0 K measured on a different crystal (crys- 
tal #2 of Fig. 9) at field H = 0. The vertical dotted red line 
indicates the T c . 



obtained on the same crystal, we also measured the low- 
T C P (T) on a different crystal #2 as shown in Fig. 11(a), 
for which the p{T) data are shown in Fig. 9(b). As shown 
in the inset of Fig. 10, a rather sharp heat capacity jump 
is observed at T c = 0.92(5) K due to the superconduct- 
ing transition, where we define T c to be the transition 
midpoint. Two important differences are observed be- 
tween the C P (T) and p(T) data on crystal #2. First, 
the T c obtained from the two measurements are differ- 
ent: T c = 0.92(5) K from C P (T) and T c = 1.5(1) K from 
p(T). Second, while the superconductivity is suppressed 
to a temperature below 0.45 K by H = 0.05 T in C P (T), 
superconductivity occurs at « 1 K even at H = 0.15 T 
in p(T). These two observations suggest the presence of 
surface sheath superconductivity in SrPd2As2 which is 
probed by p(T) whereas C P (T) measures the bulk super- 
conductivity. 

A fit of the normal-state C P (T)/T versus T 2 data using 
the H = 0.1 T data for 0.45 K < T < 4.4 K and the H = 
0dataforl.3K<T< 4.4 K in Fig. 11(b) by Eq. (3) (the 
fitted data from 3.2 K to 4.4 K are not shown) gives j n — 
6.43(3) mJ/molK 2 , j3 = 0.369(8) mJ/molK 4 and S = 
3.7(5) //J/molK 6 as shown by the red curve. The Debyc 
temperature estimated from /? using Eq. (4) is 0d = 
298(3) K. A fit of C P (T) in Fig. 10 over the entire T range 
(2-300 K) by Eqs. (5) with j n fixed to the above value 
and using the Padc approximant 54 in place of the Debyc 
function gives <dr> = 245(3) K. The fit is shown by the 
red curve in Fig. 10. The value of 0d is smaller than the 



15 



o 

E 

E^ 

c 

O 





: (a) 


1 1 1 1 1 1 1 1 1 1 








SrPd 2 As 2 


■ 


15 






A 






T =0.92K 

C 


x ■ 

y - 


10 
5 




&** H || c 


-o-0T - 

a 0.03 T " 

-v 0.05 T ! 








-n— 0.10T ■ 


n 


1 * * 


, 1 ,,,,!,, , 






0.0 



0.5 



1.0 
T(K) 



1.5 



2.0 



15 



*10 

o 
E 

E, 
H 5 

a. 

O 



1 1 1 I 

; (b) 


i i i I i i i 
SrPd 2 As 2 


| I I I | I I I 


■ « 

_ o 

°0 
" o 




^-tr 3 *^^ 


- 




H || c 




o 


Data, H = T . 


. 


D 


Data, H = 0.1 T . 


, , , i 


,,,!,,, 


-Fit 

. , 



10 



T 2 (K 2 ) 



FIG. 11: (Color online) (a) Heat capacity C p versus temper- 
ature T of a SrPd2As2 single crystal [crystal # 2 of Fig. 9] 
for 0.45 K < T < 2.0 K measured at the indicated values of 
applied magnetic field H with H \\ c. (b) C p /T vs. T 2 for 
0.45 K < T < 3.2 K with H = and 0.1 T. The red curve is 
a fit of the H = 0.1 T data for 0.45 K < T < 4.4 K and the 
H = data for 1.3 K < T < 4.4 K by Eq. (3) (the fitted data 
from 3.2 K to 4.4 K are not shown). 



value of 298(3) K obtained from the low-T fit, indicating 
a T-dependent 0d- 54 The parameters obtained from the 
analyses of the normal-state C p (T) data are summarized 
in Table IV. 

The electron-phono n coupling constant is estimated 
from Eq. (7) as A e - P h = 0.443 using /i* = 0.13, 
T c = 0.92 K and 6 D = 298 K. Then from Eq. (6) 
we obtain the bare Sommerfeld coefficient as 70 = 
4.46 mJ/molK 2 . From the value of 70 we estimate 
T>(E-p) = 1.89(1) states/ (eV f.u.) for both spin directions 
from Eq. (6), which is very close to that of CaPd2As2 
in Table IV. The Fermi velocity obtained from Eq. (9) 
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FIG. 12: (Color online) (a) The electronic contribution C e 
versus temperature T at low T obtained by subtracting the 
phonon contribution from the measured C P (T) of SrPd2As2. 
(b) C c /T versus T. The solid red curves in (a) and (b) are 
the theoretical predictions of the a-model for a = 1.30 and 
the dashed blue curve in (b) is for a = 1.50. 



is vp = 1.69 x 10 8 cm/s, and the mean free path for 
this value of vf and po = 4.9 ^tf^cm using Eq. (11) is 
I = 11.6 ran. The oj p and Al(0) estimated from Eqs. (14) 
and (16), respectively, are listed in Table VII. 

The electronic contribution C C (T) to the measured 
low-T C P (T) of SrPd 2 As 2 , obtained by subtracting (3T 3 + 
5T 5 from C P (T) according to Eq. (3), is plotted ver- 
sus T in Fig. 12(a), and C C {T)/T is plotted versus T in 
Fig. 12(b). Utilizing the entropy-conserving construction 
in Fig. 12(b) we obtain AC C (T C )/T C = 5.0(2) mJ/molK 2 , 
AC e (T c ) = 4.6(2) mJ/molK using T c = 0.92(5) K and 
AC c (T c )/7 n T c = 0.77(5) using 7n = 6.43(3) mJ/molK 2 . 
The value of AC c (T c )/j n T c is significantly smaller than 
the BCS weak-coupling value of 1.43, as was also the case 
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for CaPd.2As2 discussed above, and from Eq. (19) we ob- 
tain a = 1.30(4) which may be compared with the BCS 
value of 1.764. In Figs. 12(a) and 12(b) we show as the 
red curves the respective theoretical predictions of the a- 
model obtained using a = 1.30 in Eqs. (21). Also shown 
as the dashed blue curve in Fig. 12(b) is the theoretical 
prediction for a = 1.50, which is the value that best fits 
the H C (T) data in Fig. 13(a) below. 

The H C (T) for SrPd2As2 is obtained by integrating 
the C e (T) data in Fig. 12 according to Eqs. (22) and the 
results arc shown in Fig. 13(a). The shoulder just above 
T c arises from the high-T shoulder in the C e (T) data in 
Figs. 12(a) and 12(b). The value of H c (0) is calculated 
from Eq. (24) using a = 1.30, 7 n = 6.43 mJ/molK 2 and 
T c = 0.92 K, yielding H c (0) = 5.47 mT. The theoretical 
prediction of the a-model in Eqs. (25) for H C (T) using 
these parameters is plotted in Fig. 13(a). Although the 
T dependence of the data is reproduced, the calculated 
magnitude does not agree with the data. A better fit as 
shown in Fig. 13(a) by the blue dashed curve is obtained 
using a = 1.50, for which H c (0) is calculated as above to 
be 6.3 mT. 

The H C 2 versus T data obtained from the above 
C P (T) and p(T) measurements with H \\ c are shown in 
Fig. 13(b). As already mentioned the two measurements 
show different T c 's in zero field. Consistent with this dif- 
ference, the H C 2(T) behavior derived from the p(T) data 
indicates the presence of a surface sheath critical field 
H C 3 that is larger than the bulk H C 2. The above val- 
ues of H c (0) for both a = 1.30 or 1.50 are much smaller 
than the bulk -ff C 2(0) extrapolated from the heat capacity 
data in Fig. 13(b), indicating type-II superconductivity 
in SrPd2As2 as was also found above to be the case in 
CaPd 2 As 2 . 

From the C P (T) data in Fig. 13(b) one obtains 
dH c2 {T)/dT\ T=Tc = -0.12(4) T/K. Equation (26) 
then gives H° 2 rh \o) = 0.077 T in the clean limit and 
H° rh (0) = 0.073 T in the dirty limit. Using a = 1.30, 
the Pauli-limiting field is obtained from Eq. (28) as 
iJp(0) = 1.37 T and the Maki parameter in Eq. (29) 
is then «m = 0.08. An estimate of H C 2(T) obtained 
from the WHH prediction in Eqs. (30) using om = 0.08 
and A so = is shown as the solid red curve in Fig. 13(b), 
from which we obtain -ff C 2(0) = 70 mT. 

The Ginzburg-Landau parameter is obtained from 
Eq. (33) using H c2 (0) = 70 mT and H c (0) = 5.5 mT, 
yielding kql = 9.0. The lower critical field estimated 
from Eq. (34) is H cl (0) = 0.94 mT. The Ginzburg- 
Landau coherence length at T = obtained from Eq. (35) 
is £(0) = 69 nm which together with £ = 11.6 nm 
and Eq. (38) for the dirty limit gives £o = 509 nm. 
The Fermi velocity estimated from £o using Eq. (41) is 
t>F = 0.25 x 10 8 cm/s for a = 1.30 which is of the same or- 
der as the above estimated value of vp from the density 
of states. Corresponding values of the above parame- 
ters for a — 1.50 were also calculated. A summary of 
the measured and derived superconducting parameters 
for SrPd2As2 is given in Table VII. 
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FIG. 13: (Color online) (a) Thermodynamic critical field 
H C (T) obtained from the free energy considerations from the 
zero field superconducting state C P (T) data. The solid red 
and dashed blue curves are the theoretical predictions of the 
a-model for a = 1.30 and 1.50, respectively, (b) Upper critical 
field H C 2 (T) of SrPd2 As2 determined from the electrical resis- 
tivity p(T,H) and heat capacity C P (T,H) data in Figs. 9(b) 
and 10(b), respectively. The red curve is the prediction of the 
WHH theory in Eqs. (30) for a M = 0.07 and 0.08 and A so = 0. 



C. Magnetization and Magnetic Susceptibility 

The ZFC x(T) = M(T)/H data for a SrPd 2 As 2 single 
crystal measured in H = 3.0 T arc shown in Fig. 14 to- 
gether with the intrinsic x obtained from fitting M(H) 
isotherm data for H > 2 T in the Appendix by Eq. (A.l). 
These two data sets are in excellent agreement over the 
whole T range, indicating the near absence of ferromag- 
netic and saturable paramagnetic impurities in the crys- 
tal. The very small upturns below ~ 10 K in x{T) in 
Fig. 14 are attributed to a trace amount of paramagnetic 
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TABLE VII: Measured and derived superconducting and rel- 
evant normal state parameters for SrPd2As2. T c : bulk su- 
perconducting transition temperature; j n : observed Sommer- 
feld coefficient of the linear term in the low-T normal-state 
heat capacity; A _ p h: electron-phonon coupling constant; 70 
bare Sommerfeld coefficient; I: mean free path at low T; uj p 
plasma angular frequency; Al: London penetration depth 
A c fi: magnetic penetration depth; H®2 h '- orbital upper crit- 
ical magnetic field; A: superconducting order parameter; 
a — A(0)/kBT c ; AC C : heat capacity jump at T c ; Hp: Pauli 
limiting upper critical field; «m: Maki parameter; H c , H c i, 
H C 2'. thermodynamic, lower critical, and fitted upper critical 
magnetic fields, respectively; kgl: Ginzburg-Landau param- 
eter; £: Ginzburg-Landau coherence length; £0: BCS super- 
conducting coherence length. 
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FIG. 14: (Color online) Zero-field-cooled magnetic suscepti- 
bility x °f a SrPd2 As2 single crystal versus temperature T in a 
magnetic field H — 3.0 T applied along the c-axis (xc H || c) 
and in the afe-plane (xab, H _l_ c). The filled stars represent 
the intrinsic x obtained in the Appendix from fitting M(H) 
isotherm data by Eq. (A.l). The lines joining the stars are 
guides to the eye. 



10~ 4 cm 3 /mol. The Pauli susceptibility estimated from 
Eq. (46) using V{Ey) = 1.89 states/eV f.u. for both spin 
directions is xp = 6.1 x 10 -5 cm 3 /mol. The xl = 
—2.0 x 10 -5 cm 3 /mol was obtained by taking m* = m c in 
Eq. (47). Then (xvv) is obtained from these values using 
Eq. (45). The various \ contributions are summarized in 
Table VI. 



VI. PHYSICAL PROPERTIES OF BaPd 2 As 2 
CRYSTALS 



impurities. The \ is diamagnetic and exhibits a weak In- 
dependence with a strong anisotropy Xab > Xc- The large 
anisotropy in \ most likely originates from anisotropy in 
the paramagnetic Van Vleck orbital contribution xvv m 
Eq. (45), although the strong temperature dependence 
of the anisotropy is unusual and unexpected for such a 
compound and for which the origin is not clear. De- 
spite having the same crystal structure, the x anisotropy 
in SrPd2As2 is opposite to that of CaPd2As2 in Fig. 8 
for which we observed Xab < Xc- This difference in the 
sign of the anisotropy between the two compounds is ev- 
idently again attributable to a difference in the sign of 
the anisotropy in the Van Vleck contributions between 
the two compounds. 

The powder and temperature (1.8-300 K) average 
of the intrinsic x obtained from the M(H) isotherms 
in the Appendix is (x) = —3.9 x 10~ 5 cm 3 /mol. 
The diamagnetic core susceptibility estimated using the 
atomic diamagnetic susceptibilities 72 is Xcorc = —1-93 x 



A. Electrical Resistivity 

The in-plane p(T) data for a BaPd2As2 crystal are 
shown in Fig. 15. The T dependence of p reveals metal- 
lic behavior with a very small residual resistivity po ~ 
1 /ifJcm and a large RRR rj 25, demonstrating the high 
quality of the crystal. The scale of the resistivity for 
BaPd2As2 is smaller and the RRR is much larger than 
those of CaPd2As2 and SrPd2As2- The p(T) data were 
fitted by Eqs. (2) for 2 K < T < 300 K using the ana- 
lytic Padc approximant function 54 as shown by red curve 
in Fig. 15. The excellent fit obtained yielded the fitting 
parameters po = 1.02(1) pQ cm, p(Or) = 8.84(1) pflcm 
and Or = 114(1) K. The value of the constant TZ ob- 
tained using Eq. (2c) is 1Z = 9.34/if2cm. The fit param- 
eters are summarized in Table III. 

The expanded low-T plot of p(T) in the inset of Fig. 15 
reveals an onset of superconductivity at T consct ss 2.0 K 
for H = with zero resistance at about 0.6 K. However, 
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FIG. 15: (Color online) In-plane electrical resistivity p of a 
BaPd2As2 crystal versus temperature T measured in zero 
magnetic field H. The red curve is a fit by the Bloch- 
Griineisen model. Inset: Expanded plot of the low-T p(T) 
data showing the onset of superconductivity at T < 2.0 K 
measured at H — and 0.05 T applied along the c-axis. We 
infer that this superconductivity is filamentary and not bulk. 



FIG. 16: (Color online) Heat capacity C p of a BaPd2As2 sin- 
gle crystal versus temperature T for 2.4 K < T < 300 K 
measured in zero magnetic field. The red solid curve is the 
fitted sum of the contributions from the Debye lattice heat 
capacity CvDcbye(T) and predetermined electronic heat ca- 
pacity 7 n T according to Eq. (5a). Inset: C p /T versus T 2 for 
0.45 K < T < 5 K. The red curve is a fit of the data by 
Eq. (3) for 0.45 K < T < 5 K. 



the transition width is very large compared to those of 
CaPd2As2 and SrPd2As2, and the heat capacity measure- 
ments in the following section show no evidence for su- 
perconductivity above 0.45 K. An applied field of 0.05 T 
decreases T consot by ps 0.8 K, as shown. 



bulk C P (T) data above 0.45 K in the inset of Fig. 16, 
which indicates that there is no bulk superconductivity 
in BaPd2As2 and hence the superconductivity detected 
by the p(T) measurements is filamentary in nature. 



B. Heat Capacity 



C. Magnetization and Magnetic Susceptibility 



The C P (T) data for BaPd2As2 are shown in Fig. 16. 
Like the other two compounds, the C P (T = 300 K) = 
124 J/molK is close to the expected high-T classical 
Dulong-Petit value. An expanded low-T plot of C p /T 
versus T 2 is shown in the inset of Fig. 16. A fit of 
the data with 0.45 K < T < 5 K by Eq. (3) yields 
7„ = 4.79(2) mJ/molK 2 , /3 = 0.638(5) mJ/molK 4 and 
5 = 4.0(3) /iJ/molK 6 , as shown by the red curve in the 
inset. A V(E F ) = 2.03(1) statcs/(cV f.u.) for both spin 
directions is estimated from j n by setting A _ p h = in 
Eq. (6). A value 0d = 248(1) K is obtained from /J using 
Eq. (4). 

A value 6 D = 227(2) K is obtained by fitting the C P (T) 
data by Eqs. (5) over the entire T range (2-300 K), as 
shown by the red curve in Fig. 16. Here again we used an- 
alytic Pade approximant function 54 for CvDebyc(T) and 
set 7n to the fixed value obtained above from the fit to 
the low-T C P (T) data. The normal-state parameters ob- 
tained from these fits are summarized in Table IV. 

In contrast to the observation of the onset of super- 
conductivity in the p{T) data at w 2.0 K in the inset 
of Fig. 15, no corresponding feature is observed in the 



The anisotropic X (T) = M(T)/H data of a BaPd 2 As 2 
single crystal in H = 3.0 T are shown in Fig. 17. The 
X is diamagnctic over the whole T range and exhibits 
a weak T-depcndence. We find that Xc > Xab, which 
has the same sign of the x anisotropy as in CaPd2As2 
but is opposite to that of SrPd2As2. One must keep 
in mind in making these comparisons that the crystal 
structure of BaPd2As2 is different from that of CaPd2As2 
and SrPd2As2- The intrinsic anisotropic susceptibilities 
obtained from the slopes of high-field linear fits to the 
M(H) isotherms in the Appendix arc shown by solid stars 
in Fig. 17. These data are in rather good agreement with 
the x(T) = M(T)/H data in Fig. 17 above 25 K, and 
indicate that the low-T upturns below 25 K are due to a 
small amount of paramagnetic impurities. 

The different contributions to the intrinsic x °f 
BaPd2As2 were estimated following the same approach 
as for the Ca and Sr members above. The powder- and 
temperature-average of the intrinsic x obtained from the 
analyses of the M(H) isotherms from 1.8 to 300 K in the 
Appendix is (x) = — 1.24xl0 -4 cm 3 /mol. We also obtain 
^ corc = — 2.16xl0 _4 cm 3 /mol, xp = 6.6xl0~ 5 cm 3 /mol, 
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FIG. 17: (Color online) Zero-field-cooled magnetic suscepti- 
bility x of a BaPd2As2 single crystal as a function of temper- 
ature T in the temperature range 1.8-300 K measured in a 
magnetic field H — 3.0 T applied along the c axis (xc, H || c) 
and in the ab plane (Xab, H _L c). The filled stars are 
the intrinsic x obtained in the Appendix from fitting M(H) 
isotherm data by Eq. (A.l). The lines joining the stars are 
guides to the eye. 



and xl = -2.2 x 10" 5 
Eq. (47)]. A value (xvv, 
obtained using Eq. (45). 
marized in Table VI. 



cm 3 /mol [using m* — m c in 
= 4.8 x 10~ 5 cm 3 /mol is then 
These contributions are sum- 



VII. MAGNETIC PENETRATION DEPTH 
MEASUREMENTS OF CaPd 2 As 2 AND SrPd 2 As 2 

Figure 18 shows the temperature variation of the ab- 
plane magnetic penetration depth, AA c g (T) = \ e ff(T) — 
A c fi(0), measured in CaPd2As2 and SrPd2As2 crystals, 
represented by open circles and triangles, respectively. 
The absolute value of the penetration depth was ob- 
tained using the TDR technique by matching the fre- 
quency shift, A/(T), to the skin depth, S, calculated from 
the resistivity. The superconducting transition temper- 
ature was determined as the temperature of the maxi- 
mum of dAX g/dT. The determined T c 's are 1.34 K and 
1.26 K for CaPd2As2 and SrPd2As2, respectively. These 
values are higher than the bulk T c 's of 1.27(3) K and 
0.92(5) K determined from respective C p (T) data (Ta- 
bles V and VII, respectively). Even so, the actual onset 
of the diamagnetic response is observed at even higher 
temperatures, T° nsct = 1.50 K and 1.72 K for CaPd 2 As 2 
and SrPd2As2, respectively. The T dependences of AA e ff 
for the two compounds up to T c are shown in the inset 
of Fig. 18. 

At low temperatures, the AA ff (T) in Fig. 18 of each 
sample shows a clear saturation on cooling, which is an 
indication of a fully-gapped superconducting order pa- 
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FIG. 18: The change AA eff = A cff (T) - A cff (T -► 0) of 
the magnetic penetration depth A e ff measured in CaPd 2 As 2 
(filled circles) and SrPd 2 As 2 (filled triangles). The data for 
SrPd 2 As 2 are shifted vertically upwards by 30 nm for clarity. 
The solid curve for each compound is the best fit of the data 
by the prediction in Eq. (49) for a single-gap s-wave BCS su- 
perconductor for T < T c /3. The vertical dashed line is the 
upper temperature limit for the fits. Inset: AX e g(T) for both 
compounds up to T — T c . The data for SrPd 2 As 2 are all 
shifted upwards by 2 /im for clarity. 



rameter in both compounds. Our compounds arc dirty- 
limit superconductors, for which the magnetic penetra- 
tion depth in the single-band model for T/T c <C 1 in 
Tinkham's notation 58 is 



1 + |, 



(48a) 



where Xl(T) is the clean-limit BCS London penetra- 
tion depth prediction, yielding the T = expression in 
Eq. (37), and 



Aeff(T) = A L (T) 
Aeff(0) A L (0) " 

Defining AA(T) = X(T) - A(0), one obtains 

AA eff (T) _ AA L (T) 
Aoff(0) A L (0) 



(48b) 



(48c) 



The right-hand side of this equation is just the clean-limit 
BCS prediction for local electrodynamics given by 46,49,60 



AA L (T) 



'ttA(O) 



exp 



A(0) 



A L (0) y 2k B T 
Combining Eqs. (48c) and (49) gives 



AA cff (T) = A cff (0)i 



'ttA(O) 
2k B T 



exp 



knT 



A(0) 



(49) 



knT 



(50) 
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The experimental data are fitted well up to T sa T c /3 
by Eq. (50) as shown by the solid curves in Fig. 18, where 
the fitting parameters arc A e fj(0) = 210 ± 60 nm and 
A(0)/fc B = 2.02 ± 0.14 K for CaPd 2 As 2 and A cff (0) = 
170 ± 70 nm and A(0)/fc B = 2.05 ± 0.20 K for SrPd 2 As 2 . 
The listed errors arc systematic errors obtained from the 
spread of the fitting parameters depending on different 
choices of the upper temperature limit near T c /3. 

Using the bulk T c values in Tables V and VII and 
the above values of A(0)/A;b, we obtain a = 1.59(14) 
for CaPd 2 As 2 and a = 2.23(0.32) for SrPd 2 As 2 . The 
value of a for CaPd 2 As 2 is identical within the error bars 
to the value of 1.58(2) in Eq. (20) that was determined 
from the heat capacity jump, both of which are smaller 
than the value a B cs ~ 1.764 expected for an isotropic 
weak-coupling BCS superconductor. 46,60 This reduction 
is most likely due to a moderate anisotropy of the or- 
der parameter 46 rather than multiple order parameters, 
because well-known multi-gap superconductors such as 
MgB 2 , 75 NbSe 2 , 76 and LiFeAs (Ref. 77) have shown much 
lower values of a for the smaller gap. The accuracy of 
a for SrPd 2 As 2 that we could write down is uncertain 
because of the significantly larger superconducting tran- 
sition width obtained from the heat capacity measure- 
ments for this compound, and will therefore not be fur- 
ther considered. 

The above parameter values obtained for CaPd 2 As 2 
and SrPd 2 As 2 are listed in Tables V and VII, respec- 
tively 



TABLE VIII: Parameters obtained from fitting M(H) 
isotherms of APd 2 As 2 (A = Ba, Ca, Sr) at 1.8 K by Eqs. (A.l) 



and (A. 2 ), where 



and and Sir 



2. Here M s is 



the saturation magnetization of ferromagnetic impurities, \ 
is the intrinsic susceptibility, and /j mp is the molar fraction of 
the paramagnetic impurities. 



Compound 


field 


M s 


X 


/imp 




direction 


1 G cm 3 \ 


/ 10 -5 cm!) 


(io- 4 ) 


CaPd 2 As 2 


H 1c 


0.06(7) 


-15.38(5) 


1.01(2) 




H || c 


0.48(3) 


-14.20(3) 


0.86(1) 


SrPd 2 As 2 


H _Lc 


0.002(4) 


-2.60(2) 






H || c 


0.02(1) 


-5.91(2) 




BaPd 2 As 2 


#_Lc 


0.18(3) 


-14.59(4) 






H || c 


0.20(4) 


-11.07(3) 





tion using the same value of a. 

While the p(T) data of CaPd 2 As 2 exhibit a T c con- 
sistent with that obtained from the C P (T) data, the 
p{T) data of SrPd 2 As 2 exhibit a higher T c most likely 
due to surface sheath superconductivity. Our analysis 
of the normal and superconducting state p(T, H) and 
C P (T, H) data of these two compounds and estimated 
superconducting parameters indicate type-II supercon- 
ductivity with small thermodynamic critical fields and 
upper critical fields that are much smaller than those of 
the dopcd-FcAs based superconductors. 



VIII. 



SUMMARY AND CONCLUSIONS 



The crystallographic, electronic transport, thermal, 
magnetic and superconducting properties of APd 2 As 2 
(A = Ca, Sr, Ba) single crystals were investigated. 
The magnetic measurements in the normal state reveal 
anisotropic diamagnetism with Xc > Xab for CaPd 2 As 2 
and BaPd 2 As 2 , and \ab > Xc for SrPd 2 As 2 . The x(T), 
p(T) and C p (T) data indicate sp-band metal behavior of 
all three compounds and provide conclusive evidence for 
bulk superconductivity in CaPd 2 As 2 and SrPd 2 As 2 but 
only filamentary superconductivity in BaPd 2 As 2 which 
has a different crystal structure. The superconduct- 
ing transition temperatures as estinratcd from the zero- 
field C P (T) data are T c = 1.27(3) K for CaPd 2 As 2 and 
T c = 0.92(5) K for SrPd 2 As 2 . The heat capacity jump 
at T c , AC C (T C ), of CaPd 2 As 2 in H = is extremely 
sharp, which allows unambiguous analyses of the de- 
rived electronic contribution C C (T) in the superconduct- 
ing state. The AC c (T c )/-f n T c = 1.14(3) is significantly 
smaller than the BCS prediction of 1.43. We analyzed 
the C e (T) data in the superconducting state using the 
a-model, 46 < 61 where a = A(0)/fc B T c . A good fit to the 
data was obtained using a = 1.58, which is significantly 
smaller than the BCS value of 1.764 which we surmise 
is due to anisotropy in the superconducting s-wave gap. 
The H C (T) is also in agreement with the a-model prcdic- 
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Appendix: Presentation and Analysis of M(H) 
Isotherms 



1. CaPd 2 As 2 

The M(H) isotherms for a CaPd 2 As 2 crystal measured 
at eight temperatures between 1.8 and 300 K for H ap- 
plied both along the c-axis (M c , H \\ c) and in the ab- 
plane (M a b, H _L c) are shown in Fig. 19. Consistent 
with the X; the M is diamagnetic (negative) and exhibits 
weak anisotropy with M C (H) > M ab (H). For T > 50 K 
the M(H) curves are almost linear in H, however at low- 
T a slight nonlinearity is observed that can be attributed 
to the presence of a small amount of saturable param- 
agnetic (PM) impurities. Further, the presence of trace 
amount of ferromagnetic (FM) impurities is also inferred 
from the M(H) curves. 

We estimated the FM impurity contribution by fitting 
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FIG. 19: (Color online) Isothermal magnetization M of 
CaPd2As2 versus applied magnetic field H at different tem- 
peratures, as listed, for magnetic fields applied (a) in the 
ab-plane (H _L c) and, (b) along the c-axis (H || c). The 
solid curves are fits of the M(H) data by Eqs. (A. 2) for 
1.0 < H < 5.5 T. 



the M{H) data for T > 50 K and if > 2 T by 
M{H) = M s + X H, 



(A.l) 



where M s is the FM impurity saturation magnetization. 
For T > 50 K the M s is found to be almost T-indcpcndcnt 
and anisotropic with values for H \\ c and H _L c listed in 
Table VIII. The M s value of 0.48 Gcm 3 /mol for H \\ c 
is equivalent to the saturation magnetization of 39 mo- 
lar ppm of Fe metal impurities suggesting that the only 
trace amounts of FM impurities are present in the crystal. 
However due to the small magnitude of the diamagnetic 
X even trace amounts of FM impurities are observable in 
X and M measurements. 

Once the FM impurity contributions M s to the magne- 
tizations were determined, we analyzed the low-T M(H) 



data for both field directions by 

M(T, H) = M s + X H + ./h„ P i\4 mp Bs imp (»)• (A.2a) 

Here x is the intrinsic susceptibility of the compound, 
/imp is the molar fraction of PM impurities, M Sitn = 
A'A3impA f B>S'i m p is the PM impurity saturation magnetiza- 
tion, Na is Avogadro's number, ^b is the Bohr magneton, 
and 3i mp and ^imp are the spectroscopic splitting factor 
(p-factor) and the spin of the impurities, respectively. 
Our unconventional definition of the Brillouin function 



D 



Sin 



IS 



TS 



Bs lmp (x) 



2Si, 



(2S imp + l)coth (25 imp + l) 



where 



— coth 



^B (T — Simp 




2 J 
(A.2b) 

(A.2c) 



A Weiss temperature #i mp is included in the argument of 
Bs im (x) to take into account for interactions between the 
paramagnetic impurities in an average mean-field way. 
While fitting the M(H) data the impurity (/-factor was 



set to <7i, 



2 and the M s values for H _L c and H 



were set to the values listed in Table VIII. The M(H) 
data for both H directions were fitted for magnetic fields 
in the range 1.0 < H < 5.5 T. The S[ mp values for both 
H directions were found to be Si mp = 2.0(2); therefore, 



in the final fits we set S'imp = 2. 



The #i mp values for 



both field directions were found to be close to zero and 



therefore in the final fits we set 9 n 



0. The solid 



curves in Fig. 19 show the final fits of the M(H) data 
by Eqs. (A. 2). The parameters obtained from the fits 
of the M(H) isotherms at T = 1.8 K are listed in Ta- 
ble VIII. The intrinsic x values obtained from the fits 
of the M(H) data at different temperatures are shown 
by stars in Fig. 8. The T dependence of the intrinsic 
X values clearly indicates that the low-T upturns in the 
measured x(T) = M(T)/H data are extrinsic. 



2. SrPd 2 As 2 

The isothermal M(H) data for a SrPd2As2 crystal 
at different T are shown in Fig. 20. Similar to the 
X(T) = M{T)/H data in Fig. 14, the M(H) data ex- 
hibit anisotropic diamagnetic behavior with M a b(H) > 
M C (H). In order to obtain the contributions from the 
FM impurities the M(H) data for if > 2 T at T > 1.8 K 
were fitted by Eq. (A.l) which gave the average M s values 
of 0.002(4) for H _L c and 0.02(1) for H || c. The fitting 
parameters for both H \\ c and H _L c for T = 1.8 K 
are listed in Table VIII. The intrinsic susceptibilities 
obtained from the analysis of M(H) data by Eq. (A.l) 
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FIG. 20: (Color online) Isothermal magnetization M of 
SrPd2As2 versus magnetic field H at different temperatures 
as in Fig. 19. 



FIG. 21: (Color online) Isothermal magnetization M(H) of 
BaPd2As2 versus magnetic field measured at different tem- 
peratures as in Fig. 19. 



are shown by stars in Fig. 14. The paramagnetic im- 
purity concentration is zero within our resolution, since 
no curvature in the M(H) curves on cooling to low T in 
addition to that due to the FM impurities was detected. 



metal impurities, respectively. As in the Sr compound, 
the paramagnetic impurity concentration is zero within 
our resolution, since no curvature in the M(H) curves 
on cooling to low T in addition to that due to the FM 
impurities was detected. 



BaPd 2 As 2 



The M(H) isotherms for a BaPd2As2 crtstal at differ- 
ent T are shown in Fig. 21. The M(H) curves exhibit 
weakly anisotropic diamagnetic behavior with M C (H) > 
M a b(H). The intrinsic x was obtained by fitting the 
M(H) isotherms at each T by Eq. (A.l) for H > 2 T 
which yielded a temperature-averaged (for T > 50 K) 
FM saturation values Mf = 0.18(3) G cm 3 /mol and M s c 
= 0.20(4) Gcm 3 /mol which arc equivalent to the mag- 
netization contributions from 12 and 16 molar ppm of Fc 
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